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INTRODUCTION 

The  diffusion  equation  of  Physics  has  been  used  to  analyze  unsteady  heat 
transfer,  boundary  layer  velocity  distribution,  long  line  electrical  voltage 
fluctuation,  and  salt-solute  penetration.  The  general  mathematical  expression 
is 

3T/3t  =  a  32T/3X2  , 

where  particular  physical  constraints  determine  the  context.  There  are  two 
approaches  to  the  problem  statement  and  solution;  the  one  most  widely  used 
being  transformation,  and  the  other,  finite  difference  techniques,  employ 
variations  of  summing  averages  of  term  values  established  by  unique  methods. 
This  report  considers  an  averaging  type  solution  in  algebraic  format.  The 
final  result  consists  of  a  series  of  discrete  polynomials  with  rational 
coefficients  which  describe  the  dependent  variable  at  each  time-distance 
coordinate  in  the  manner  of  the  non-reflecting  Schmidt  plot. 

PROCEDURE 

Essentially,  the  differential  equation  is  recast  as  a  finite  difference 
expression  which  is  transposed  first  to  geometric  and  then  to  polynomial 
algebraic  form.  The  polynomials,  representing  discrete  solutions  to  the 
differential  equation,  are  analyzed  by  differencing  techniques  whereby  the 
numerical  coefficients  of  common  diagonal  terms  are  found  to  be  expressible  in 
a  generalized  matrix. 

From  the  one  dimensional  partial  differential  equation 


8T 

3t 


(1) 


where  T  is  the  dependent  variable,  t  is  an 
independent  variable,  and  a  is  a  constant, 
temperature,  t  is  time,  A  is  distance,  and 
of  the  finite  difference  procedure  gives 


independent  variable,  K  is  an 
In  heat  transfer  problems,  T  is 
a  is  the  diffusivity.  Application 


a  At  .2 


(2) 


with  At  the  finite  difference  in  time,  Ax  the  finite  difference  in  distance, 
A  T  the  time  variable  effecting  a  change  in  T,  and  A  T  the  distance  variable 
effecting  a  change  in  T.  * 

By  expansion  the  equation  becomes 


T  (t  +  1) 
n _ 

At 


(t) 


(n+1 ) 


-  2  T 


(t) 


+  T 


(t) 


(n-1) 


(3) 


L 
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where  subscript  n  refers  to  the  x  increments  and  superscript  t  refers  to  the  t 
increments.  Schmidt1  developed  the  graphical  form  shown  on  Figure  1  with  the 
stepwise  linear  temperature  gradients  across  adjacent  layers  of  material. 

Since  the  change  in  internal  energy  within  a  layer  of  material  over  a  finite 
time  is  the  difference  between  the  heat  flow  in  and  heat  flow  out,  the  ^ 
corresponding  temperature  increment  becomes  a  function  of  the  ratio  &x 
and  it  is  convenient  to  select  this  ratio  as  unity,  leading  to  2aAt 


Also,  a  geometric  simplification  results  from  defining  the  graphical 
proportions  as 

.  «> 

where  the  Axq  is  the  graphical  analogue  of  the  initial  film  coefficient. 

Table  1  shows  the  discrete  algebraic  expressions  for  the  time-temperature- 
distance  intersections  of  Figure  1 .  Along  the  diagonals  of  Table  1  a  matched 
power  polynomial  appears  and  the  coordinate  expression  for  T  in  time  and 
distance  takes  the  general  form 


T  (N,P)  = 


(A  -  Bm  +  Cm 


±  mk) 


where  N  is  a  distance  index,  P  is  a  time  index, 


(P  +  N)  -  2  - 

2 


sin  (P  +  N)  tt/2 


(the  external  denominator  exponent),  k  is  the  terminal  exponent  of  m,  j  is  k 
minus  the  exponent  of  m  (the  individual  term  denominator  exponent),  and 

v  =  m  (T0°  -  Tj0) 

with  A,B,  C  ...  numerical  coefficients  of  the  interior  terms  of  the  equation. 

To  establish  T  (N,  P)  for  any  $  and  m  (which  include  the  physical  constraints) 
it  is  necessary  to  establish  the  precise  values  of  the  coefficients  A,  B,  C 
and  this  is  the  object  of  the  current  investigation. 

Any  full  expression,  T  (N,P),  can  be  developed  from  a  Gregory-Newton 
formulation^  of  the  separate  terms  A,  B,  C,  for  the  bounded  diffusion  equation 
as  shown  in  Figure  2,  Some  interesting  progressions  do  result  but  an  alternate 
and  more  geometric  presentation  is  available  from  the  direct  finite  difference 
tables . 


2  Sutton,  G.P.,  Rocket  Propulsion  Elements,  John  Wiley  &  Sons,  New  York  1956. 

Spiegel,  M.R.,  Theory  and  Problems  of  Finite  Differences  and  Finite 
Difference  Equations ,  Schaum's  Outline  Series  in  Mathematics,  McGraw-Hill  Book 
Company,  New  York,  etc.,  1971,  pp.  36-44. 


Table  2  is  also  extracted,  by  difference  equation  procedure,  from  Table  1 
and  generates  the  coefficients  of  the  constant  term,  A,  for  the  coordinate 
expression  of  distance  and  time.  The  starting  point  is  within  the  heavy  box 
of  column  7,  These  numbers,  17548,  25147,  35401,  49024  and  66868,  were  found 
by  direct  calculation  using  Figure  1  and  Table  1  and  are  the  constant 
numerator  terms  only.  Table  3  lists  the  complete  polynomial  expressions  for 
these  coordinates.  By  the  usual  differencing,  the  7th  through  the  1st  columns 
are  established.  It  is  then  possible  to  work  vertically  using  the  regression 
in  column  2,  back  to  zero;  and  then  to  complete  the  elements  of  columns  1 
through  6.  Noting  the  resulting  bias  progression  at  the  tops  of  these 
columns,  the  next  step  is  continue  diagonally  (I,  II,  III)  to  column  8  and, 
using  column  7  as  a  summation,  verify  the  vertical  sequence  of  column  8. 
Columns  9,  10,  11,  ...,  are  generated  similarly. 

Within  the  individual  frames  containing  the  coefficients  is  a 
parenthesized  pair  of  numbers  which  indicate  the  distance  and  time  coordinate. 
These  indices  run  diagonally  upwards  at  constant  distance  and  bi-sequentially 
as  time.  Tables  4  through  9  are  formed  by  the  same  procedure  and  extend 
arbitrarily  to  the  6th  power  of  m.  However,  a  different  sequence  appears 
along  diagonals  I,  II  and  III  according  to  the  power  of  m.  Table  10 
summarizes  this  behavior  and  reveals  yet  another  correlation,  shown  mainly  bv 
column  4,  from  which  the  adjacent  columns  can  be  constructed  adinfinitum.  The 
final  coincidence  occurs  from  a  reinspection  of  Tables  2,  and  4-9  where  the 
digital  vertical  counting  column  (1,  2,  3,  4,  5  ...)  conjoins  the  power  of  tn 
and  the  time  sequence  of  the  first  distance  diagonal  (1,3),  (1,5),  (1,7), 
etc.,  by  an  interval  of  3  in  the  counter  according  to  Table  11. 

RESULTS 

To  write  any  term  defining  the  dependent  variable  in  time  and  distance, 
Tables  2  through  11  are  used  to  form  the  numerical  coefficients  in  Equation 
(6), 


T  (N,P)  =  - —  f  A  -  Bm  +  Cm2  -  .  .  .  +  .  •  .  ±  mK  ] 

2  2J 

For  T  (1,15)  for  instance, 

N  =  1, 

P  =  15, 


k 


P  -  N  -  1  sin  (P  -  N)  tt  /2  | 

2 


(15-1)-  j  sin  7  it ] 
2 


P  +  N  -  2  |  sin  (P  +  N)  tt/2  1 


(15  +  1)  -  2 
2 


«  7 ,  and 


j  =  (k  -  term  exponent  of  m)  *  (7  -  terra  exponent  of  m) 
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The  finite  difference  tables  for  the  (1,15),  j  variables  are  then 
reconstructed  (Tables  2,  and  4  -9)  using  Tables  10  and  11  and  the  respective 
numerators  determined,  yielding 


whereby 


Term 

Numerator 

Exponent 

»i  j*» 

Value 

of  V 

A 

51480 

0 

7 

B 

39796 

1 

6 

C 

20264 

2 

5 

D 

7050 

3 

4 

E 

1672 

4 

3 

F 

260 

5 

2 

G 

24 

6 

1 

H 

1 

7 

0 

hence 


T  0,15) 


= 

128 


51480 

128 

260m5  + 


39796m  20264m2 

64  32 


7050m3  1672m4 

16  8 
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TEMPERATURE 
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1.  Schmidt  Diagram 


TABLE  1.  TIME  LOCATION 


N 

P  - — 

1 

2 

3 

4 

3 

6 

i 

i 

0 

0 

(0/2) (3  -  m) 

(0/2) (3  -  m) 

(0/a) (30/a  -  9/2  m  +  m2) 

(0/4)  (30/4  -  9/2  m  +  m2.: 

2 

0/2 

0/2 

(0/a) (7/2  -  m) 

(0/4 )  (7/2  -  m) 

(0/8) (38/4  -  10/2  m  +  m2 

3 

0/A 

0/a 

(0/8)(a  -  m) 

(0/8) (4  -  m) 

4 

'Sl 

CD 

0/8 

(0/16) (9/2  -  m) 

5 

0/16 

0/16 

6 

0/32 

7 

8 

1L 


TABLE  1.  (continued) 


_ 7 _ 

(0/8) (35/2  -  29/2  m  +  12/2  <n2  -  m3) 
(0/8) (35/4  -  10/2  m  +  m2) 
(0/16)(47/4  -  11/2  m  +  m2) 

(0/16) (9/2  -  m) 

(0/32) (5  -  m) 

0/32 

0/64 


_ 8 _ 

(0/8) (35/2  -  29/2  m  +  m2  -  m3) 

(0/16)  (187/8  -  69/4  m  +  13/2  tn2  -  m3) 
(0/16) (47/4  -  11/2  m  +  m2) 

(0/32) (57/4  -  12/2  m  +  m2) 

(0/32) (5  -  m) 

(0/64) (11/2  -  m) 

0/64 

0/128 
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TABLE  1.  (continued) 


_  9 _ 

(0/16) (630/16  -  325/8  m  +  95/4  m2  -  15/2  m3  +  m*) 
(0/16)(187/8  -  69/4  m  +  13/2  m2  -  m3) 

(0/32) (244/8  -  81/4  m  +  14/2  m2  -  m5) 

(0/32) (57/4  _  12/2  m  +  m2) 

(0/64) (68/4  -  13/2  m  +  m2) 

( 0/64 ) (11/2  -  m) 

(0/128) (6  -  m) 

0/128 

0/256 
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TABLE  1.  (continued) 


_ 10 _ 

(0/16) (630/16  -  325/8  m  +  95/4  m2  -  15/2  m3  +  m4) 
(0/32) (874/16  -  406/8  m  +  109/4  m2  -  10/2  m3  +  m4) 
(0/32) (244/8  -  81/4  m  +  14/2  m2  -  m3) 

(0/64) (312/8  -  94/4  m  +  15/2  m2  -  m3) 

(0/64) (68/4  -  13/2  m  +  m2) 

(0/128) (80/4  -  14/2  m  +  m2) 

(0/128) (6  -  m) 

(0/256) (13/2  -  m) 

0/256 

0/512 


17 


TABLE  1.  (continued) 


_ 11 _ 

(0/32) (1386/16  -  843/8  m  +  312/4  m2  -  141/4  m3  +  18/2  m4  -  m5) 
(0/32) (874/16  -  406/8  m  +  109/4  m2  -  16/2  m3  +  m4) 

(0/64) (1186/16  -  300/8  m  +  124/4  m2  -  17/2  m3  +  m4) 

(0/64) (312/8  -  94/4  m  +  15/2  m2  -  m3) 

(0/128) (392/8  -  298/4  m  +  16/2  m2  -  m3) 

(0/128) (60/4  -  14/2  m  +  m2) 

(0/256) (93/4  -  15/2  m  +  m2) 

(0/256) (13/2  -  m) 

(0/512) (14/2  -  m) 

0/512 

0/1024 


18 


TABLE  1.  (continued) 


_ 12 _ 

(0/32) (1386/16  -  843/8  m  +  312/4  m2  -  141/4  m3  +  18/2  m*  -  m5) 

(0/64) (1979/16  -  1093/8  m  +  374/4  m2  -  79/2  m3  +  19/2  m*  -  m5) 

(0/64) (1186/16  -  500/8  m  +  124/4  m2  -  17/2  m3  +  m*) 
(0/128)(1578/16  -  608/8  m  +  140/4  m2  -  18/2  m3  +  m*) 

(0/128) (392/8  -  298/4  m  +  16/2  m2  -  m3) 

(0/256) (485/8  -  123/4  m  +  17/2  m2  -  m3) 

(0/256) (93/4  -  15/2  m  +  m2) 

(0/512) (214/8  -  16/2  m  +  m2) 

(0/512) (14/2  -  m) 

(0/1024) (15/2  -  m) 


0/1024 


TABLE  1.  (continued) 


_ 13 _ 

(0/64) (3003/16  -  4165/16  m  +  1841/8  m2  -  532/4  m3  -  196/4  m4  +  21/2  m5  -  m6) 
(0/64) (1979/16  -  1093/8  m  +  374/4  m2  -  79/2  m3  +  19/2  m4  -  m5) 

(0/128) (2762/16  -  1397/8  m  +  444/4  m2  -  88/2  m3  +  20/2  m4  -  m5) 

(0/128) (1578/16  -608/8  m  +  140/4  m2  -  18/2  m3  +  m4) 

(0/256) (2063/16  -  731/8  m  +  157/4  m2  -  19/2  m3  +  m4) 

(0/256) (485/8  -  123/4  m  +  17/2  m2  -  m3) 

(0/512) (592/8  -  139/4  m  +  17/2  m2  -  m3) 

(0/512) (214/8  -16/2  m  +  m2; 

(0/1024 ) (244/8  -  17/2  m  +  m2) 

(0/1024) (15/2  -  m) 

(0/2048) (16/2  -  m) 

0/2048 

0/4096 
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TABLE  1.  (continued) 


_ 14 _ 

(0/64) (3003/16  -  4163/16  m  +  1841/8  m2  _  532/4  m3  +  196/4  m4  -  21/2  m5  +  m6) 
(0/128) (4387/16  -  5562/16  m  +  2285/8  m2  -  310/2  m3  +  216/4  m4  -  22/2  m5  +  m6) 

(0/128) (2786/16  -  2794/16  m  +  444/4  m2  -  88/2  m3  +  20/2  m4  -  m5) 

(0/256) (7599/32  -  3525/16  m  +  1045/8  m2  -  195/4  m3  +  21/2  m4  -  m5) 

(0/256) (2063/16  -  731/8  m  +  157/4  m2  -  19/2  m3  +  m4) 

(0/512) (2655/16  -  870/8  m  +  175/4  m2  -  20/2  m3  +  m4) 

(0/512) (592/8  -  139/4  m  +  18/2  m2  -  m3) 

(0/1024) (714/8  -  156/4  m  +  19/2  m2  -  m3) 

(0/1024) ( 244/8  -  17/2  m  +  m2; 

(0/2048) (276/8  -  18/2  m  +  m2) 

(0/2048) (16/2  -  m; 

(0/4096) (17/2  -  nr 
0/4096 
0/8192 


TABLE  1.  (continued) 


_ 15  _ 

(0/128) (6435/16  -  9949/16  m  +  5066/8  m2  -  3525/8  m3  +  836/4  m4  -  130/2  m5  +  24/2  m6  -  m7) 

(0/128) (4387/16  -  5562/16  m  +  2285/8  m2  -  310/2  m3  +  216/4  m4  -  22/2  m5  +  m6) 

(0/256) (25147/64  -  14649/32  m  +  5615/16  m2  -  1485/8  m3  +  237/4  m4  -  23/2  m5  +  m6) 

(0/256) (7599/32  -  3525/16  m  +  1045/8  m2  -  195/4  m3  +  21/2  m4  -  m5) 

(0/512) (10254/32  -  4395/16  m  +  1220/8  m2  -  215/4  m3  +  22/2  m4  -  m5) 

(0/512) (2655/16  -  873/8  m  +  175/4  m2  -  20/2  m3  +  m4) 

(0/1024) (3369/16  -  513/4  m  +  97/2  m2  -  21/2  m3  +  m4) 

(0/1024) (714/8  -  156/4  m  +  19/2  m2  -  m3) 

(0/2048) (852/8  -  174/4  m  +  20/2  m2  -  m3) 

(0/2048) (138/4  -  18/2  m  +  m2) 

(0/4096) (155/4  -  19/2  m  +  m2} 

(0/4096) (17/2  -  m) 

(0,  8192  (18/2  -  ra. 

0/8192 

0/16384 


22 


TABLE  !•  (continued) 


r 


_ 16 _ 

(0/128) (6435/16  -  9949/16  m  +  5066/8  m2  -  3525/8  m3  +  836/4  m4  -  130/2  m5  +  24/2  m6  -  m7) 
(0/256) (76627/128  -  43947/64  m  +  25879/32  m2  -  8485/16  m3  +  1090/8  m4  -  283/4  m5  +  23/2m6-m 
(0/512) (35401/64  -  9502/16  m  +  6835/16  m2  -  825/4  m3  +  259/4  m4  -  24/2  m5  +  m6) 

(0/512) (10254/32  -  4395/16  m  +  1220/8  m2  -  215/4  m3  +  22/2  m4  -  m5) 

(0/1024) (13623/32  -  5421/16  m  +  1414/8  m2  -  236/4  m3  +  23/8  m4  -  m5) 

(0/1024) (3369/16  -  513/4  m  +  97/2  m2  -  21/3  m3  +  m4) 

(0/2048) (4221/16  -  600/4  m  +  107/2  m2  -  22/2  m3  +  m4) 

( 0/2048 )( 852/8  -  174/4  m  +  20/2  m2  -  m3) 

(0/4096) (1007/8  -  193/4  m  +  21/2  m2  -  m3) 

(0/4096; (155/4  -  19/2  m  +  m2, 

(0/8192) (173/4  -  20/2  m  +  m2, 

(0/8192 . (18/2  -  m 
( 0/16384 ) ; 19  2  -  m 
0/16384 
0/32768 


23 


21 


(•reyory-Mowton  Trans pos it  i 
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TABLE  10.  LINE  PROGRESSION  CORRELATIONS 
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TABLE  11.  VERTICAL  COUNTER  CORRELATION 


Term 

Exponent 

"IV"  Column 
Value 

Distance-Time 

Coordinate 

0 

6 

(1,3) 

1 

9 

(1,5) 

2 

12 

(1,7) 

3 

13 

(1,9) 

4 

18 

(1,11) 

3 

21 

(1,13) 

6 

24 

(1,15) 

7 

27 

(1,17) 

8 

30 

(1,19) 

9 

33 

(1,21) 

9 


ackncwlkdgemeht 


The  finite  difference  and  difference  equation  path  followed  in  this 
report  must  be  ascribed  to  Mr,  Stowell  Elder,  who  also  graciously  expended 
time,  effort  and  encouragement  to  complete  the  development. 
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LIST  or  SYMBOLS 


a  constant  (dif fusivity ) 

h  exponent  of  2  in  the  external  denominator 

j  exponent  of  2  in  each  term  denominator 

k  exponent  of  ,rmM  in  final  term 

m  constant  (— - ^ -  ) 

Ax  *  Axq 

n  number  of  x  increments 

t  independent  variable  (time) 

X  independent  variable  (distance) 

numerical  coefficients 
N  distance  index 

P  time  index 

T  Dependent  variable  (temperature) 

3  indicating  partial  d i f f erentr at ive 

A  indicating  difference 

A^T  time  variable  effecting  a  change  in  T 

AT  distance  variable  effecting  a  change  in  T 

<;  external  numerator 
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